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Abstract 

We consider T-shaped, two-dimensional quantum waveguides containing attractive or repulsive 
impurities with a smooth, realistic shape, and study how the resonance behavior of the total 
conductance depends upon the strength of the defect potential and the geometry of the device. 
The resonance parameters are determined locating the relevant 5-matrix poles in the Riemann 
energy surface. The total scattering operator is obtained from the S'-matrices of the various 
constituent segments of the device through the -^"-product composition rule. This allows for a 
numerically stable evaluation of the scattering matrix and of the resonance parameters. 

PACS numbers: 73.63.Nm Quantum wires, 73. 23. Ad Ballistic transport, 72.10.Fk Scattering by point 
defects, dislocations, surfaces, and other imperfections (including Kondo effect) 



Preprint DFPD/07/TH11 



1 



I. INTRODUCTION 



The discovery of conductance quantization in microcostrictions has prompted a 

great deal of theoretical and experimental activity on electron transport in such systems 
3|, 14, l5|. It has been soon realized that the presence of an impurity in the waveguide can 
modify the shape of the conductance dramatically, especially near the thresholds where 
propagation modes are opened. These effects have been investigated in many model calcu- 



ations, for both point-like defects 



a lid, 



11 



a, 13, |9| and more realistic, finite-range defect shapes 



12j. The scattering of the electron wave off the impurity produces resonance 



phenomena in the conductance; in particular, for attractive impurities one has deep trans- 
mission minima, which can be interpreted as due to the coupling of a propagation mode in 
one subbband with a localized state formed from another subband. This resonant suppres- 
sion of ballistic conductance is similar to the Fano resonances one observes in atomic and 
nuclear physics, and has been the subject of detailed numerical 1(3, 11, 12] and theoretical 
HQ investigations. Needless to say, whereas theoretical arguments can give insight into 
the general mechanisms underlying the formation of resonances, only detailed numerical 
experiments can provide information about the dependence of resonance parameters upon 
the strength, shape, and position of the impurity potential inside the quantum waveguide. 

The step-wise structure of conductance can be deeply influenced also when the waveguide 
boundary is modified, with respect to an idealized, straight shape. In particular, the effects 



of a sidearm, or stub, attached to the duct have been the subject of severa 



3, 



16. 



13, OS, 09, 



20 



2l| since the pioneering papers of Sols et al. 



22 



investigations 



231 ]. where a 



T-shaped waveguide was considered as a possible candidate for transistor action based on 
quantum interference. The conductance as a function of the Fermi energy exhibits deep 
transmission minima, which can be again interpreted as reflection resonances due to the 



excitation of quasi-bound states in the cavity region [19] . A great deal of theoretical activity 
has been devoted to disentangle the physical factors influencing Fano line shapes in these 
mesoscopic systems. Indeed, Fano resonances are interesting in themselves, because they are 
extremely sensitive to the details of the scatte ring process, and may be exploited to probe 
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271 ] : at the same time, their strong 



the degree of coherence in the scattering device 
dependence upon the energy of the scattered particle make them promising candidates for 
practical applications, such as the spin filtering action when spin degeneracy is removed by 
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the application of an external magnetic field or by use of the Rashba effect 28j]. At low 
level densities, and for systems coupled to very few channels, quantum interference between 
quantum states plays a prominent role, and a detailed quantum-mechanical study of the 
scattering device is required. The dependence of the spectral properties of quantum dots 
attached to quantum waveguides upon the dot's dimensions has been the subject of detailed 
analysis in recent years [29I, Q . 

To the best of our knowledge, up to now an analysis of the simultaneous effects of im- 
purities and stubs in quantum interference devices has never been accomplished. It is the 
purpose of this paper to provide a first contribution, in order to fill this gap. In particular, 
we have studied how the resonance features of the conductance depend upon the strength 
of an embedded defect, and the size of the stub. To this end, we have determined how the 
resonance poles of the scattering operator for the device move on the multi-sheeted Riemann 
energy surface, as the device parameters are varied. For stubbed waveguides with realistic 
defects, this requires an accurate numerical solution of the two-dimensional Schrodinger 
equation, in a large range of values for the stub's dimensions. This is made possible com- 
bining mode-matching techniques with an S"-matrix approach. The device is regarded as a 
cascade of uniform waveguide sections, connected through interfaces or junctions, where a 
change in the transverse dimensions occurs, or the electron experiences a variation in the 
strength of the interaction. The scattering parameters of the device can be then determined 

juilding blocks through the -^-product 
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321 • Since forward and backward 



starting from the scattering properties of the various 
composition rule for the partial scattering operators 
states are treated separately by the S'-matrix, so that the propagating and less localized 
components play a prominent role, numerical stability is guaranteed as some dimension 
of the system becomes large. At the same time, a varying number of basis functions can 
be accommodated in a natural way in the scattering matrix approach; this is particularly 
important in non-uniform waveguides, where one expects that a different number of wave- 
function components has to be considered as the transverse size of the device changes. Thus, 
the employment of the scattering operator, together with the ^--product composition rule, 
opens the way to an accurate description of devices with a fairly large number of differ- 
ent constituent segments, as required when dealing with impurities having a smooth profile 
along the propagation direction. They are in fact approximated by a cascade of several thin 
layers, in which the dependence upon the propagation coordinate is neglected. 
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The paper is organized as follows. In Sect. [IT] we review the scattering matrix approach 
to ballistic transport in quantum waveguides. In particular, we point out the compatibility 
between the -^■-product composition rule and the non-trivial block-wise structure of the 
partial S-matrices, when a different number of basis functions is chosen in different slices 
of the device. With reference to this point, we would like to point out that, whereas the 



numerical stability of the S'-matrix has been widely recognized in the literature [3|, |3l|], its 



flexibility with respect to the number of wave function components has been less stressed 



in the past, with the exception of Ref. 33|. In Sect. II III we give several examples of how 
the device parameters influence the pole location in the energy-plane, and how this affects 
the resonant behavior of the conductance. Our main conclusions are briefly summarized in 
Sect. HVl 



II. SCATTERING FORMALISM FOR QUANTUM TRANSPORT IN STUBBED 
WAVEGUIDES 

We consider the T-shaped device illustrated in Fig. [TJ It consists of a uniform guide 
of indefinite length and width b, with a sidearm (stub) having width c and length l s . We 
shall assume simple hard wall boundary conditions, so that the two-dimensional electron 
wave function has to vanish on the boundary of the device. The stub may contain a region, 
where there is a defect, or an external applied field. This region, to which we shall simply 
refer as "defect", is exhibited by the shaded area in Fig. [TJ The electron's wave function 
y) is then the solution of the two-dimensional Schrodinger equation for a given total 
energy E, the electron in the conduction band being endowed with an effective mass m*. 
In the present paper we have chosen for m* the value 0.067m e , which is appropriate for the 
Al^Gai-^As/GaAs interface. The potential field due to the defect will be represented by an 
interaction term V(x,y), whereas the effect of the confining potential is taken into account 
through the boundary conditions that y) vanishes along the edge of the straight duct 
and of the stub. Because of the presence of the potential term V(x,y), the full Schrodinger 
equation is not in general separable. As is well known [5], one can reduce it to an (in 
principle) infinite set of one-dimensional differential equations by expanding its solution 
into the transverse mode eigenfunctions (fin (y) in the leads (r = I) and in the stub (r = s). 
Under the assumption of hard-wall boundary conditions, these basis functions are just the 
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FIG. 1: A stubbed quantum waveguide of width b and infinite length, with a stub of width c and 
length l s . The stub contains a defect with dimensions w x l^. 

eigenfunctions of infinite square wells in the transverse direction, with widths b and c in the 
leads and in the stub, respectively. 

The whole device can be divided into five regions; the left and right ducts, the two 
empty cavities in the sidearm, and the defect. The total wave function in the different 
regions of the stub can be piecewise expanded in terms of the basis functions <f)n\ Once 
this expansion is inserted into the Schrodinger equation, one immediately finds that the 
expansion coefficients satisfy uncoupled differential equations, and can be written in terms 
of forward and backward propagating waves in the regions of the stub where the potential 
V(x, y) is negligible; the corresponding wave numbers k$ are related to the total energy E 
and to the transverse eigenenergies by 



Us) 



2m* /mr\ 2 



2m* (E - e ( n 

w 



(2.1) 



These relations essentially express the conservation of energy; once the n-th transverse 
eigenmode has been excited, only the energy (h 2 /2m*)kn^ is left for the electron propagation 
along the x-direction. If E < €n\ the wave number has to be taken as purely imaginary, 
i.e., kn^ = iKn\ with 

/ /T7-7T \ 2 9m* 

(2.2) 



K 



(s) 



fnn\ 2 2m* 



The associated propagating waves become then real, exponentially decaying or growing wave 
functions, and one has an evanescent mode, or closed channel. 

In presence of the defect a somewhat more involved procedure is required, since one gets 
coupled equations for the expansion coefficients. These equations can be straightforwardly 
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uncoupled, and the coefficients written in terms of forward and backward propagating waves, 
if one assumes for the potential V(x, y) the simple form 

V(x,y) = V f(y)e(x-x 2 )e(x 3 -x) , (2.3) 

where Q(x — Xi) is the Heaviside step function. In such a case one arrives at the required 



result diagonalizing the Hamiltonian 



h 2 d 2 

2m* dy 2 



for the transverse motion, in the 



model space spanned by the basis functions (f>n\y) [6|. In the potential region of the stub the 
expansion coefficients can then be written as linear superpositions of plane-wave components 
e ±ik { n \x-x 2 ) w ith wa ve vectors 



*j» = ^WEM , (2 .4) 

the quantities £ n being the eigenvalues of the transverse Hamiltonian. Here also, if E < £ n , 
the wave number is purely imaginary, in completely analogy to Eq. (|2.2p . and one has 
evanescent, non propagating modes in presence of the defect interaction. 

Similar considerations can be repeated for the wave function in the leads, the only dif- 
ference being that now the transverse eigenmodes have to be used. The wave vectors 
in the leads are given by expressions quite similar to Eqs. (12.11) and (12.21) for open and 

(s) 

closed channels, respectively, with the stub transverse eigenenergies e« replaced by the lead 
corresponding quantities e>n = - — ( — ] . For the time being, we shall not specify the 



2m* \ b 

actual boundary conditions the wave function has to satisfy asymptotically, and we let the 
presence of incoming and outgoing waves from both the left and from the right in the ducts. 
We just limit ourselves to observe that, in correspondence to closed channels, one must ex- 
clude incoming waves from the left or from the right, in order to avoid divergent components 
in the wave function. 

The unknown coefficients appearing in the expansion of the wave function in the various 
regions can be related to one another by matching the wave function and its first derivative 
at the various interfaces. The linear transformation yielding the coefficients on one side 
of an interface in terms of those on the other side represents what is usually referred to 
as the transfer matrix It is well-known that the transfer matrix suffers from several 



limitations 
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341 ] . First of all, its generalization to the case where there is a different 
number Nl of basis functions on the left of the matching surface with respect to the right 
is not straightforward. This is particularly relevant at the lead/stub or stub/lead interfaces, 



where a different number of open and closed channels may occur in the cavity with respect 
to the lead. A second problem arises in presence of evanescent modes, whose occurrence 
makes the inversion of the transfer matrix numerically troublesome when one evaluates the 
transmission coefficients as the dimensions of the device get large. The above shortcomings 
can be all overcome resorting to the scattering matrix 3|, |4|. In such a case, one expresses 
the outgoing amplitudes in terms of the incoming ones, to get 
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(2.5) 



where the sub-matrices are given by 



1.1 

Sff = F (A 21 - BaaBr/Au) , (2.6a) 

Sg = F(B 22 + B 21 Br 1 1 B 12 ) , (2.6b) 

S$ = F (A n + A 12 AJ 2 1 A 21 ) , (2.6c) 

Sg = F (-B 12 + A 12 A^B 22 ) , (2.6d) 



with 



F= (A 22 + B 21 B U 1 A 12 ) 1 (2.7a) 
F= (Bh + Ai 2 A 2 ~ 2 1 B 2 i)^ 1 . (2.7b) 

The actual structure of the matrices Ay and B^- depends upon the matching conditions 
one is considering. Whatsoever these conditions may be, one does not need to assume the 
same number of basis functions on the two sides of the matching surface, since the building 
blocks of the total S'-matrix are not required to be square. As a matter of fact, if Nl wave 
function components are retained on the left, whereas Nr basis functions are considered on 
the right, the whole sub- matrix Sn has dimensions Nl x Nl, and S 22 is a square Nr x Nr 
array, whereas S12 and S 2 i have dimensions Nl, x Nr and Nr x Nl, respectively. Overall, 
the total S'-matrix expresses the Nl + Nr outgoing amplitudes in terms of the Nl + Nr 
coefficients of the incoming waves. As for the diagonal matrix T+\ it takes into account 
the propagation of the wave function in the slice on the left of the matching surface, and is 
defined according to 

\ A + fnm — C "nm , V^ -0 / 



7 



where Ax is the width of the slice, and k^ are the corresponding propagation wave numbers. 
It is worth to observe that in the scattering matrix only the phase factors e +lk " )Ax appear. In 
correspondence to an evanescent mode, one actually has the exponentially decaying factors 
e Kn x , which are no longer mixed up with the growing functions e Kn x . As the numerical 
calculations confirm, these features make the whole formalism numerically stable even when 
some characteristic length of the system, and in particular the stub's width c, gets large. 

Finally, we give the actual structure of the matrices Ay and By at the lead/stub and 
cavity/defect interfaces. For the former one has An = A 12 = J, B n = B 12 = 1, and 
A 2 i = A 2 2 = K^, B 2 i = B 2 2 = JK( S \ The N R x N L matrix J contains the overlap integrals 
among the different basis functions on the left and on the right of the discontinuity, 

* 4 ] (y)<t>$(y)dy , (2.9) 

J is the transpose of J, while and are Nl X Nl and Nr x Nr diagonal matrices, 
whose diagonal elements represent the wave vectors k$ and kn , in the lead and in the stub, 
respectively. For the cavity/defect interface, on the other hand, one has A n = A 12 = 1, 
Bn = B 12 = Z, and A 2 i = A 22 = K^, B 2i = B 22 = ZK< d >, where Z is the matrix 
diagonalizing the transverse Hamiltonian in the interaction region, and the diagonal matrix 
j£( d ) contains the wave vectors kn defined by Eq. (12.41) . For the sake of simplicity, we have 
assumed an equal number of basis functions on the left and on the right of the cavity/defect 
boundary; as a matter of fact, the potential V(x, y) in general is not so strong, as to perturb 
the transverse eigenvalues in an essential way, and good convergence can be achieved with 
an equal number of wave function components outside and inside the defect. 

The unperturbed propagation along the slice preceding the interface can be put in better 
evidence by re- writing Eq. (12.51) as a suitable combination of e. 
This can be achieved through the -^--product of S"-matrices 



ementary scattering operators. 



3|, |3ll, |32| , which expresses the 



overall scattering matrix S in terms of the partial scattering matrices and as 

S S 

S = | 11 12 ] = S (a) *S (b) , (2.10) 

S21 S22, 



s 



where 





°11 °22 




22 a ll 



(a) Q (6) 




(2.11b) 



(2.11d) 



(2.11a) 



(2.11c) 



It is worth to observe that this composition rule does not require the sub-matrices S^- 
and to be square matrices. In particular, the S-matrix appearing in Eq. (12.51) can be 
factorized as follows 



where the block-wise anti-diagonal factor takes into account the free propagation 
through the slice on the left of the interface, and the reduced S-matrix exhibits the 
effects on the wave function of the matching at the interface between the two considered 
regions of the device. The former is a 2Nl x 2Nl matrix, whereas the latter has dimensions 
(N L + N R ) x (N L + N R ), as detailed above. 

Eqs. (I2.10p and (12.121) have been systematically employed to evaluate the total S-matrix, 
starting from the partial S-matrices at the lead/stub and cavity/defect boundaries, as well 
as from the scattering operators for the free propagation of the electron along the various 
slices of the device. As noted above, while there are physical reasons to choose a different 
number Ni of transverse basis functions in the lead with respect to the sidearm, one can 
choose the same number N s of components in the empty slices of the stub and inside the 
defect. Thus, one has an (Ni + N s ) x (JVj + N s ) partial S-matrix S(/|s) at the lead/stub 
boundary, and a (2N S x 2N S ) partial S-matrix S(c|d) for the cavity/defect propagation. One 
can verify by inspection that, when the two operators are combined through the composition 
rule ( 12.101) . one gets again a matrix with dimensions (Ni + N s ) x (jVj + N s ). Things are a bit 
different when one considers the stub/lead interface on the right of the device. We shall not 
enter here into the detailed structure of the stub/lead interface S'-matrix. We limit ourselves 
to note that, once the lead/stub S'-matrix S(Z|s) < ' / ' > has been evaluated, the corresponding S'- 
matrix S(s|/) ( ' / ' ) on the right can be obtained without further computational effort, through 




(2.12) 
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the correspondence rule 

/s( s |/) (7) 
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u) 



(s(l\s)$ 8{l\s)S\ 



\S(J|s)g> S(l\s)[\>J 



u) 



(2.13) 



A similar result holds for the S'-matrices at the cavity/defect and defect /cavity interfaces. 
As a consequence of Eq. (12. 13ft . once the (JVj + N s ) x (A^ + N s ) partial S'-matrix describing 
the propagation inside the device is combined with the (N s + N t ) x (N s + Ni) S-matrix S(s\l) 
through the -^--product composition rule (I2.10p . one gets an overall 2Ni x 2Ni scattering 
operator S^, which effectively takes into account the flux percolation through open and 
closed channels in the sidearm. 

Potentials having a non-trivial x dependence can be accommodated in our approach 
through a slicing technique 12J, ll7|]. This is particularly easy to do for interactions which 
are factorisable in their x and y dependence, such as the double Gaussian 

V(x,y) = v e- p2{x - Xc)2 - a2{v - Vc)2 . (2.14) 

If xq and xf are two points, where the value of the potential is negligible, we divide the 
interval [xq, xp\ into iV sub-intervals of equal width Ax = (xp — xq) /N, and mimic the 
smooth potential V(x, y) by a sequence of pseudo-defects, each of width Ax and strength 

V (i) = V e-P 2{x *- Xc)2 ( 



0. 



AO 



with Xi = Xq + iAx. The conductance of the waveguide in presence of the smooth potential 
(12. 14=11 is then evaluated for increasing N, until convergence of the results is obtained. 

Let us finally fix the actual boundary conditions under which we are solving the scattering 
problem. We shall assume that an electron impinges on the device moving from the left, 
with no flux coming from the right. If one assumes an incoming wave of unit flux in a given 
propagation mode, then the matrix elements [Syp ) represent the reflection coefficients 
toward the left from the initial channel m into the final propagation mode n, while the 
quantities ( S^p ) give the transmission amplitudes to the right from mode m into mode 



n. Finally, for a two-probe device the total conductance G can be evaluated through the 



Biittiker formula 



2e 



.0) 



n — hjn 



AT) 
5 21 



(2.15) 



where the summation is restricted to the open channels in the leads. 



10 



III. SELECTED NUMERICAL EXAMPLES 



Our aim is to study the dependence of resonance parameters in a stubbed quantum 
wave guide upon the strength of an embedded defect and the geometry of the device, with 
particular reference to the stub height. We shall do this by locating the relevant S-matrix 
poles in the multi-sheeted energy plane, with cuts starting along the real axis from the various 
scattering thresholds. From now on, all lengths are measured in units of the waveguide width 

(l\ H 2 /7T \ ^ 

b, and energies in terms of the waveguide's fundamental mode e\ = ( — ) , so that 

2m* \ b / 

the various thresholds are E n = 1, 4, 9, . . .. Sheets will be specified giving the sign of the 
imaginary part of the lead's momenta in each channel [36']. Thus, referring to a four-channel 
situation, the symbol (+ + ++) denotes the physical sheet, where the imaginary parts of 
all the channel momenta are positive, whereas on sheet ( — h ++) one has Im/c^ < 0, 
and Imkf^ > for the other three channels. Our calculations refer to a duct with a double 
Gaussian defect described by Eq. f |2. 14|) . In all cases we have chosen x = x /b = 0, and 
Xf — 1, so that the region where the potential acts has a length 1^ = 1. Moreover, we have 
set w = 0.3; for a wave guide of width b = 100 A , the interaction region is therefore 30 A 
wide and 100 A long. The defect is placed at Y = 0.1 from the lower edge of the guide. 
The decay constants along the transverse and the propagation direction have been fixed at 
a = ba = 15 and /3 = 10, respectively, so as to ensure that the potential is negligible outside 
the interaction region. We have checked the convergence of the slicing method, and found 
that quite stable results are obtained with iV = 10 -j- 15 sub-intervals. 

For a waveguide attached to a stub of height c = c/b and length l s = l s /b = 1, one can 
expect that a larger number of basis functions is required in the stub than in the guide to 
achieve convergence, this difference increasing with increasing c. As a matter of fact, for a 
given electron's energy E there are less open channels in the waveguide than in the stub, 
the thresholds in the two regions being Ejp = n 2 and E n s ' = En /c 2 , respectively. It follows 
that, for a given energy E, the number of open channels in the stub scales as c\^E, whereas 
in the guide it scales as \^E. Since a sensible calculation has to include at least all the 
propagating modes, a larger number of basis functions is required in the stubbed region. 
As for the evanescent modes, one may recall that a closed-channel component is the less 
effective the shorter its attenuation length £ n = 1/ (bn n ) 17] . The attenuation lengths in 
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the lead and in the cavity are given in terms of the total energy E by 

= I i £n = / = • (3-1) 

7iVn 2 -E n J^y_ E 

From Eqs. (13.11) . one has that two evanescent modes n and m, in the external region and 
in the stub, respectively, have the same attenuation length if the corresponding quantum 
numbers are related by m = cn. As a consequence, if 3 -j- 4 components are required in 
the lead, ~ 12 closed channels have to be included in the stub, to assure convergence. Our 
numerical results are consistent with these estimates for both the conductance and the pole 
locations in the complex energy plane. We have obtained stable results with four components 
in the lead, and 10 -j- 12 basis functions in the stub for 1 < c < 4. In these conditions, the 
position of the resonance poles in the .E-plane can be guaranteed with an accuracy of the 
order 10~ 5 . Only the bound-state poles turned out to be more sensitive to the number of 
employed basis functions in the interaction region, which was to be expected. In this latter 
case, we employed four channels in the duct and up to 20 modes in the internal region. 

The conductance of stubbed waveguides is characterized by a more or less complex pattern 
of transmission maxima and minima when considered as a function of the electron's energy. 
In Fig. [2] we give the conductance of a waveguide having c = 2, with a repulsive defect of 
strength V = 20, in the energy region corresponding to the first two subbands (1 < E < 9). 
One notes two transmission zeros in the first subband, and a non-trivial oscillatory behavior 
with two large dips in the second one. We looked for poles of the S'-matrix in the unphysical 

sheets ( — h ++) and ( h+) and found two resonance poles at E ~ 1.75 — 0.19i and 

E ~ 3.31 — 0.26i on the former sheet, and two pairs of neighboring poles in the ( h+) 

sheet. The value of their real parts is exhibited by the vertical dashed lines in Fig. [2j The 
two poles close to the first subband, together with the two transmission zeros, concur to 
produce the strong variations in the conductance one observes there. As for the poles on 

the ( h+) sheet, the first pair corresponds to the two almost overlapping peaks one 

sees just above the lower edge of the second subband, whereas the other two poles give 
rise to the broad resonances at the top of this energy region. It is worth to note that the 
maxima of the conductance are not aligned with the pole positions, which are very close 
to the transmission zeros in the first subband. The conductance behavior in the lowest 
subband is consistent with what one observes in simple, single-channel models for stubbed 
waveguides, where transmission resonances arise because of the presence of zero-pole pairs 

12 




FIG. 2: Conductance (in units 2e 2 /h) of a waveguide with a stub c = 2 high and a double-Gaussian 
defect of strength V = 20, in the energy region 1 < E < 9. The vertical lines show the value of 
the real part of the pole positions in the complex energy plane. The poles lie in the ( — h ++) and 
( h+) sheet in correspondence to the first and second subband,respectively. 



in the complex energy plane 371 ] . In the higher subband, one clearly has a strong coupling 

between the occurring resonances. 

We have studied how the resonance poles move as the defect strength varies with respect 

to the unperturbed situation, with the interaction switched off. Limiting ourselves to the 

first subband, for the sake of simplicity, we give the outcome of our calculations in Fig. [3], 

for a stub of height c = 2 and the strength varying in the interval —60 < V Q < 200. For each 

resonance pole in the fourth quadrant of the energy plane, one observes a complex conjugate 

pole in the first quadrant. This is a confirmation of the correctness of our calculations. 

Indeed, the multichannel S"-matrix, when re gard ed as a function of the channel momenta 

t I 

ki, has the basic property S(ki) = S*(—k*) [38], which implies that to a pole located at 
E = E&> there corresponds a pole at the complex conjugate energy E = E&*. As the defect 
strength becomes more and more repulsive, both resonance poles in the fourth quadrant move 
upwards in energy, starting from the empty-stub positions, marked by the full dots in Fig. [3] 
At the same time they move away from the real axis; the presence of a repulsive interaction 
displaces the resonance peaks at higher energy, as it was to be expected. The broadening 
of the resonances with increasing defect strength is consistent with the calculations of Ref. 
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2.5 



3.5 



Re E 



FIG. 3: Trajectories of the 5-matrix poles in the ( — h ++) sheet of the complex energy plane 
for —60 < V < 200 and a stub's height c = 2. The full and empty dots mark the empty-stub 
positions of the resonance poles and their complex conjugates, respectively. 



25| , where the effects of an impurity near the opening of a quantum dot was considered, in 



an attempt to simulate the non-monotonic resonance width found in a previous experiment 



241 ] . This effect can be attributed to the deflection of more and more electrons into the dot, 



as V increases, so that the coupling between the waveguide and the cavity gets stronger 



251 ]. Stated in equivalent terms, one may say that the wave function trapping in the region 



where the waveguide and the duct match decreases at the defect becomes more repulsive. 
In the range of strengths we have explored (0 < V < 200), on the other hand, we found 
no evidence of a non-monotonic behavior of the resonance width. For V > 200 the pole at 
higher energy moves above the upper edge of the first subband E 2 , and becomes a shadow 
pole As for the low-lying pole, we have pushed the strength V up to 1000, finding 

no signal of a decreasing width. The reason for this difference with respect to [25( can be 
traced back to the different role played by backscattering in the two cases. In Ref. [25j] the 
impurity is near the opening of the dot, whereas in our case it is embedded in the lower 
half of the waveguide. The pole trajectories in the complex energy plane depend in a rather 
sensitive way upon the position of the defect in the device, as we have found in previous 
model calculations [40:]. If the defect is displaced closer to the dot, backscattering can have 
an increased role, and a non-monotonic behavior of the resonance width may emerge. 
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For negative strengths, the two poles become closer and closer to the real axis; the first 
pole touches the real axis at E ~ 1.45 for V ~ —15.6, while the higher pole trajectory 
touches it at E ~ 2.76 when V ~ —8.4. Within the accuracy of the numerical calculations, 
we found that each pole lies just on the real axis in these conditions, at energies in the 
continuum of the first transmission band. This is an indication that for V ~ —8.4 and 
V ~ —15.6 one has a bound state embedded in the continuum (BIC) in the first subband. 
For increasingly negative strengths, the S-matrix poles move away again from the real 
energy axis. The behavior of the low-lying pole is particularly worth of attention. As the 
defect becomes more and more attractive, its trajectory bends towards the real axis, which 
is reached for V ~ —39.6, at E ~ 0.93. At this point, the complex conjugate pole coming 
from the first quadrant collides with the resonance pole, giving rise thereby to a double pole 
below the scattering threshold E n — 1. What is happening can be most clearly perceived in 
the complex plane, as exemplified in Fig. HI As the absolute value of the strength V D 

increases, the pole in the fourth quadrant of the momentum plane moves downwards towards 

(i) * 

the imaginary axis, the corresponding pole at — k\ in the third quadrant doing the same; 
the two poles collide at \mk^ ~ —0.83, when V Q ~ —39.6. Decreasing the strength further, 
one has two anti-bound-state poles, one moving downwards along the negative imaginary 
axis, the other moving upwards, until, for V ~ —46.5 it crosses the real axis, and passes 
into the upper half of the complex /c-plane as a bound state pole. Had we looked at the 
corresponding poles in the energy plane, we would have seen two poles moving on the real 
axis on the unphysical ( — sheet, one of the poles reaching the first scattering threshold 
Ei, and going back as a bound-state pole on the physical (+ + ++) sheet. This situation is 
strongly reminiscent of what happens in standard potential scattering theory, and has been 
found also in idealized, one-dimensional models of stubbed waveguides 40j. 

We have studied how the bound-state energy changes when the strength of the defect 
or the stub's height c varies. The possibility of bound-state solutions for an empty stub 
has been known since several years ago {5]. For an empty stub, the increasing of c has a 
binding effect, the binding energy Ebs becoming substantially independent upon the stub's 
size for c > 2. This can be easily explained, since the bound-state wave function is mainly 
concentrated in the region where the duct matches the stub, and for large cavities a change 
in the transverse dimension cannot affect the wave function in a substantial way. When 
the interaction is switched on, one observes that the bound state appears for increasing 
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FIG. 4: Motion of a pair of bound-state/anti-bound-state poles in the complex /cj^-plane with 
varying strength in a stubbed waveguide with c = 2. The arrows display the movement direction 
as the potential becomes more and more attractive. 

values of c for repulsive defects; for a given c, the binding energy is shifted upwards. In 
correspondence to an attractive defect, the (E BS , c) curve is shifted downwards, and a 
bound state is possible even for a smooth wave guide, when the potential is able to support 
a bound state by itself. 

The asymmetric line shapes one observes in the conductance of quantum waveguides are 
generally parametrized in terms of the Fano function 



T P = T ( 



1 + e 2 



where To is the amplitude of the Fano resonance, e = yE — E^j /f represents the reduced 
energy, and the pole location E^' in the energy plane has been written as E^ = E^ — iT. 
The Fano parameter q is a measure of the ratio between the resonant and non-resonant 
transmission amplitudes, and is related to the asymmetry of the line shape. For strictly 
one-dimensional systems, q can be evaluated starting from the positions E^°' and E^> of 
the S'-matrix zeros and poles in the energy plane, i.e, 41] 

E(R) _ £(o) 

q = ? . (3.2) 

Eq. 13.21 can be used in the present multi-channel situation below the second scattering 
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FIG. 5: Fano parameter q for the two resonances in the first subband for < V < 100 in a 
waveguide attached to a stub with c = 2. 

threshold, since in the first subband only one propagation mode really contributes to the 
conductance. To this end we have studied how the zeros of the scattering amplitude move 
as the strength varies. As expected, as the interaction becomes more and more repulsive the 
zeros move towards the upper threshold, whereas they are displaced towards lower energies 
for more and more negative strengths. In Fig. [5]we plot the resulting q parameter for the two 
resonances observed in the first subband, when the strength varies in the range < V < 100 
and c = 2. The Fano parameter of the low-lying resonance is always positive, whereas for 
the second resonance it turns out to be negative. This implies that the two line shapes 
have opposite asymmetries for all the considered strengths, as confirmed by the calculated 
conductance. In both cases the absolute value of the asymmetry parameter remains always 
different from zero, with a rather smooth dependence upon V in a rather large range of 
strength values. 

We have also studied how the position of the S'-matrix poles and zeros depends upon 
the stub's height c. The "binding" effect of an increasing c observed for the bound-state 
energies is clearly perceivable in the scattering region also. As the stub becomes longer, 
poles and zeros move towards the lower edge of the first subband, the trajectories of the 
first and second pole moving away from the real axis for c < 1.5 and c < 2.2, respectively. 
As c is increased further, the two poles come closer and closer to the real axis, giving rise to 
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FIG. 6: Fano parameter q (upper panel) and width V (lower panel) for the first (solid line) and 
second (dashed line) resonance in the first subband for 1 < c < 5 and V Q = 20. 



narrower and narrower resonances, until they leave the scattering region. We looked for the 
possible presence of BICs with varying c. No positive-energy states with a vanishing width 
emerged from our calculations, both in presence of an impurity and for an empty stub. BICs 

'or various values of the total transverse width have been found on the other hand in Ref. 

3o| , for a quantum billiard with a transverse hyperbolic profile coupled symmetrically to a 
quantum waveguide. There, arguments have been also given in favor of a rather widespread 
presence of BICs in mesoscopic systems. The fact that such zero-width states do not appear 
in the present calculations in a rather large range of stub's heights c (1 < c < 5) can be 
explained by the geometry we have chosen. Our stub represents a rather narrow resonator, 
with l s = 1; the quantum billiard considered in Ref. 30], on the other hand, has l s = 
4 and c ~ 4 -4- 5. The possible presence and location of BICs is strictly related to the 
cavity's eigenenergies, whose values and spacings depend in a crucial way upon the cavity's 
dimensions. As a matter of fact, we verified that for an empty stub symmetrically coupled to 
the leads, and of dimensions comparable to the quantum billiard of Ref. 30J there is a BIC 
at E ~ 1.420 for c = 4.6, in remarkable agreement with [30], given the different shape of the 
confining potential. We checked the effect of a host impurity on the presence and energy of 
BICs. To this end, we introduced a double-Gaussian defect centered in the cavity, having 
Id = w = 1. We found that BICs are a rather "robust" feature of the system, since for 
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Vo = 20 the zero-width state survives the switching on the interaction. When the strength 
is increased up to V = 100, one has a BIC at a slightly higher energy (E BIC — 1.465), 
provided that the stub's height is adjusted to c = 4.55, which means a variation of the 
order of 1%. Bound states in the continuum may appear also in more complicated systems. 
Recently, they have been found for two open quantum dots, coupled through a connecting 
bridge 42j. In this case also we verified that when the strength of a host interaction is 
increased from up to 100, the BIC survives, being slightly displaced upwards in energy, at 
price of a small adjustment in the neck's length l b . More precisely, one has E BIC — 3.30, 
lb = 3.59 for V = 0, and E B ic — 3.36, It, = 3.73 for V = 100. Both for the single and the 
coupled quantum dots the underlying mechanism, as viewed in the complex energy plane, is 
the same; when the strength of the interaction changes, both the pole and the transmission 
zeros move away from their original positions and are shifted upwards in energy; since they 
move at different speeds, however, if they coincided for a given value of the strength, they 
do not coincide any longer for the new strength. Their relative position can be adjusted 
to zero by a change of some characteristic length of the system, so as to have a BIC at a 
somewhat higher energy. 

In analogy with Fig. [51 in Fig. [6] we plot the Fano parameter q (upper panel) and the 
width T (lower panel) for the considered poles as functions of c. Note that the q parameter 
for the higher-energy resonance is not given for c < 1.6; as a matter of fact, for stubs 
of smaller width, this zero-pole pair cannot develop any longer in the considered energy 
region. For stubs having c < 1.6 one actually finds that the considered S"-matrix pole in the 
( — h ++) sheet is located above the scattering threshold E 2 , thereby becoming a shadow 
pole. The asymmetry parameters turns out to be rather sensitive functions of c, and can 
change sign as c varies in between 1 and 5. As a consequence, one expects a rather marked 
variation in the dependence of conductance upon the energy for stubs of different width, an 
expectation which is confirmed by our calculations. The origin of Fano lineshape reversal 
has been studied in a model coupled-channel calculation in Ref. j^j]. There, the sign change 
of q for the resonances appearing in the electron transmission spectrum of a quantum dot 
has been attributed to the coupling among shape and Feshbach resonances originating from 
open and closed channels; this effect turned out to be very sensitive to slight modifications 
of the quantum-dot geometry. A similar mechanism is presumably at work here, since as c 
increases, the poles move in the energy plane and new propagating modes open up in the 
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cavity region; as a consequence, with varying c, the q parameter of each resonance changes 
sign, and the Fano profiles one observes in the first subband may or may not have the same 
asymmetry, depending upon the stub's height one is considering. 

Finally, we observe that the analysis of the pole trajectories can be extended to the 

( h+) sheet, in correspondence to the second subband. Similarly to what has been 

found below the second threshold, one finds a strong repulsive effect with decreasing c. 
Poles, determining a resonance behavior for a given c, move up in energy until they cross 
the third threshold and become shadow poles; at the same time new poles may come into 
play, moving from lower energies. 



IV. CONCLUSIONS 



We have studied the poles in the complex energy plane of the multi-channel S'-matrix 
for a waveguide attached to a resonant cavity. We have allowed for the presence of a defect 
in the waveguide. Our main aim was to to get a quantitative insight into the dependence 
of the resonance parameters upon the defect strength and stub's height. We found that a 
decrease of the stub's size as an increase of the positive strength has a repulsive effect on the 
poles, which move in both cases to higher energies in the .E-plane. In particular, as the stub 
gets shorter, the S'-matrix poles pass the upper edge of the considered subband, becoming 
thereby shadow poles; at the same time, poles at lower energies come up, cross the lower 
threshold of the subband, and are therefore able to influence the resonant behavior of the 
conductance. For attractive defects, poles in the first and fourth quadrant move towards 
the real the strength is made more negative; with reference to the lowest energy 

region, there are indications that bound states embedded in the continuum occur for critical 
values of the strength. As the defect becomes more and more attractive, resonance and 
anti-resonance poles collide on the real axis, giving rise to anti-bound state poles moving 
on opposite directions, until a bound state pole appears on the physical sheet, in close 
analogy with what happens in one-dimensional models of stubbed waveguides [4^]. For 
repulsive defects, zero-width scattering states appear for both single and coupled quantum 
dots, provided that the resonating cavity has the proper geometry, and represent a "robust" 
phenomenon with respect to the presence of an embedded interaction. 

For the lowest subband, where only one open channel contributes to the conductance, one 
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can relate the parameters appearing in the Fano function for the line shape to the poles and 
zeros of the scattering amplitude in the energy plane. The asymmetry parameter q turns 
out to be a rather smooth function of the strength, whereas it depends more sensitively 
upon the stub's size, and may even change sign with varying c. As a consequence, one may 
or may not observe a reversal of the Fano line shape for the resonances in the first subband, 
depending upon the value of the stub's height. 

Some issues deserve certainly further investigations. The role of shadow to dominant pole 
transitions in explaining the changes of the total conductance near thresholds as c varies 
has been discussed elsewhere |44|. There is evidence of strong resonance overlapping in 
higher subbands. A quantitative study of resonance coupling effects, and their dependence 
upon the dynamical and geometrical parameters of the quantum device can be performed 
resorting to modern effective-Hamiltonian techniques [2] , in analogy to what has been done 
for microwave billiards [46|. These issues are currently under consideration. 
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